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We show how to differentiate the minimal supersymmetric extension of the standard model from 
black hole events at the Large Hadron Collider. Black holes are simulated with the CATFISH gen- 
erator. Supersymmetry simulations use a combination of PYTHIA and ISAJET. Our study, based on 
event shape variables, visible and missing momenta, and analysis of dilepton events, demonstrates 
that supersymmetry and black hole events at the LHC can be easily discriminated. 
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I. INTRODUCTION 



■ At CERN's Large Hadron Collider (LHC) fl protons will soon collide at an astonishing 800 million times per 
\ second to provide experimental evidence for the Higgs @, ||, supersymmetry (SUSY) 0, Q or extra dimensions 
[H0H1- SUSY is widely considered to be one of the best candidates for physics beyond the Standard Model (SM). It 
provides an explanation for the Higgs mass problem, a candidate for cold dark matter, and unification of low energy 
gauge couplings by introducing superpartners to SM fields (see Ref. [9( and references therein). An alternative to 
SUSY is given by phenomenological extra-dimensional models such as Large Extra Dimensions (LEDs) warped 
braneworlds [7( or universal extra dimensions [8]. Scenarios with LEDs are specially appealing. In these models, 
gravity becomes strong at the TeV scale, where radiative stability is achieved. The fundamental scale of gravity, 
M* ~ 1 TeV, is related to the observed Planck scale, Mj>\, by the relation Mpj ~ V n M™ +2 , where V n is the volume of 
the extra n-dimensional space. One of the most astounding consequences of the existence of extra-dimensions would 
be the production of subatomic Black Holes (BHs) in particle colliders [l(J, [ll| and cosmic ray showers [l2| ■ (For 
reviews, see Refs. [l3|.) 

The ATLAS [l4| and CMS [l5| experiments at the LHC are entrusted with the task of studying events with 
large transverse momentum (Pt), a signature common to both SUSY and extra dimensions. While we wait for 
these experiments to start collecting data, it is worthwhile to look into means of distinguishing SUSY and extra- 
dimensional models [l6j . Comparisons of SUSY and universal extra dimensions/little Higgs models in colliders have 
been investigated by various authors [TtJ • Discrimination of SUSY and BH events by means of dilepton events was 
recently discussed by the authors in Ref. [Hj]. In this paper, we revisit that analysis and extend it to include event 
shape variables, missin g tr ansverse momentum P T and visible energy BH and SUSY events are simulated with the 
BH generator CATFISH [19[ and the high-energy event generator PYTHIA [2fJ, respectively. SUSY masses are set with 
ISAJET [2l[. The analysis below will show that SUSY and BH events can be clearly distinguished at the LHC. BH 
events tend to be more spherical than SUSY events because of the isotropic nature of BH decay. Thus event shape 
variables, such as sphericity, provide good discriminators. On the contrary, visible energy and P T are less effective 
discriminators because of the presence of invisible channels in both SUSY and BH models, which make the amount 
of P T comparable in the two scenarios. The dilepton invariant mass is also an excellent discriminator; the SUSY 
invariant mass shows a sharp cutoff at ~ 100 GeV, which is absent in the BH model because most of the dileptons 
originate from uncorrelated events. 

The remainder of this paper is organized as follows. In Sect. [LTland lllll we briefly review the fundamentals of SUSY 
and TeV BHs which are needed for our analysis, respectively. Simulations are described in Sect. HVl The analysis of 
visible/missing momentum and event shape variables is presented in Sect. fVAl and the discrimination of SUSY and 
BH using dileptons is discussed in Sect. IV 51 Conclusions are presented in Sect. I VII 
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II. SUPERSYMMETRY AT THE LHC 

The Minimal Supersymmetric extension of the Standard Model (MSSM) Q is the simplest SUSY model. According 
to the MSSM, all SM fermions (bosons) must have a bosonic (fermionic) partner. Superpartners have identical masses, 
charges and quantum numbers of their SM counterparts, differing only in their spin. The MSSM allows for the 
unification of electromagnetic, weak and strong forces at Mgut ~ 10 16 GeV. Since we do not observe superpartners 
of SM particles at low energies, SUSY must be a broken symmetry. The SUSY breaking scale, i.e. the mass scale at 
which we expect the first SUSY particles to appear, is generally assumed to be around 1 TeV. A method of SUSY 
breaking which is mediated by gravitational interactions is supergravity (SUGRA). In its minimal version, mSUGRA 
is determined by a point in the five-dimensional moduli space with parameters: 

• mo, the common scalar mass at Mgut] 

• ™i/2j the common gaugino mass at Mgut', 

• Aq, the common trilinear coupling at Mgut] 

• tan /3, the ratio of the vacuum expectation values of the two Higgs fields; 

• /i, the sign of the Higgsino mass parameter. 

mSUGRA parameters for five typical LHC points are given in Table U 22]. Neutralino (x°), gluino (g) and squark (q) 
masses are determined by too and mi/ 2 as Xi ~ TO i/2/2, X2 ~ Xi ~ TO i/2; 9 ~ ^ m i/2 an d m(q) ~ (mo 2 + Gm 2 ^) 1 / 2 

TABLE I: Parameters for the five mSUGRA points discussed in the text. The scalar mass and the gaugino mass are given in 
GeV. 
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Visible energy, missing transverse momentum and sphericity for the five LHC points of Table [J are shown in Fig. [Tl 
where all SUSY processes except SM Higgs production have been implemented. Sparticle production at point D (open 
blue circles) is higher as squarks and gluinos are lighter. This point is usually taken as the comparison point between 
the LHC and other experiments, e.g. Tevatron [2J| and NLC (2f|. For the purposes of our analysis, the difference 
between the five points is not significant and any of them can be chosen as SUSY benchmark. In the following, 
we will consider point A. This is justified by the fact that point A allows for SUSY Higgs production [26[. Since 
BHs may evaporate into Higgs (see Sect. IIIII below ). a meaningful comparison of SUSY and BH events requires the 
presence of the Higgs channel in both models. Moreover, distinguishability of SUSY and BH events must be assessed 
by minimizing the differences between the two models. Since BH events are characterized by up to several TeV of 
missing transverse momentum, SUSY points with large such as point A, must be considered. 

A symmetry of the MSSM is i?-parity 

where B (L) is the baryon (lepton) number and s is the particle spin. All SM particles have Pr = +1 whereas their 
superpartners have Pr = —1. i?-parity implies that SUSY particles are always pair produced from SM particles. 
If i?-parity is conserved, the endpoint of a SUSY process at the LHC is a state with SM particles and two lightest 
stable SUSY particles (LSPs), which are generally neutralinos. Being colorless and chargeless, the LSPs escape the 
detector and are the source of missing transverse momentum, a leading signature of SUSY events. If i?-parity is not 
conserved, the missing transverse energy is reduced by the LSP decay. In the following, we will assume that i?-parity 
is conserved, in agreement with the MSSM (mSUGRA) scenario. 

We end this section with a list of dominant SUSY interactions at LHC point A and the definition of invariant mass. 
This is important for the following analysis because it enables us to select processes that could serve as potential 
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FIG. 1: Comparison of visible energy (left), missing transverse momentum f T (middle) and sphericity (right) for 10000 events 
for the five LHC points of Table|T](A: black filled circles, B: red filled squares, C: green filled triangles, D: blue open circles and 
E: pink open squares). 




FIG. 2: The top three SUSY decay chains and their branching ratios. 



discriminators. The third decay chain in Fig. [2] is specially interesting because it allows the separation of isolated 
leptons from the hadronic background (27j . The invariant mass is defined as 



Mi 2 = V(£i+£ 2 ) 2 -( Pl +p 2 ) 2 = V2PiP 2 (l-cos( 



where 9 is the angle between the two particles. The method of constructing invariant masses from SUSY decay chains 
has been traditionally used to calculate sparticle masses [26j . In Sect. IVBI we will use the invariant mass of isolated 
dilcptons with large Pt as a potential SUSY/BH discriminator. (For a review on lepton production at colliders, see 
Ref. [23].) In the rest frame of the second lightest neutralino, the dilepton invariant mass is 
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Since the momentum of the LSP is not constrained, the invariant mass distribution shows an edge at ~ 100 GeV. 



III. BLACK HOLES AT THE LHC 



In LED scenarios, pp collisions at the LHC could produce TeV-mass BHs with characteristic lifetimes of 10 
seconds [ljj]. Numerous studies have focused on BH signatures at the LHC [ll| and various Monte Carlo generators 
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are available for simulation purposes I19l. 123 . |30 | 1 . A quick look at BH production at the LHC reveals the following. 
According to Thome's hoop conjecture [33], a BH of mass M is formed when an object is compressed in all directions 
such that 

C < 2 7T R S (M), 

where C is the circumference of the region where the object is compacted into and R s is the Schwarzschild radius for 
a BH of mass M. An upper limit on the BH mass is obtained by assuming no gravitational energy loss at formation, 
corresponding to the black disk (BD) cross-section <jbd = vri?^. A more realistic model assumes that all the CM 
energy is not available for BH formation, some being lost as gravitational radiation (see Ref. [331 ] for a more detailed 



discussion). To estimate the energy loss, the colliding particles are treated as two Aichelburg-Sexl shock waves 34 1; 
the overlap of the shock waves forms a trapped-surface (TS) which sets a lower limit to the mass of the BH 35 1. 
(For an alternative estimate of the collisional gravitational loss, see Ref. [36].) The cross-section at the LHC involves 
summing up the contributions from all the initial partons. The cross section in the TS scenario is 

(Tpp^BH(s,n) = / 2zdz I dx j — fi(x' ,Q)fj(x/x' \Q) F a BD (xs,n) , 

ij 







where Q is four-momentum transfer squared, fi(x, Q) are the parton distribution functions, z is the normalized impact 
parameter, and F is a form factor. The cutoff in x is related to the minimum allowed mass of the object, M m i n , 
and the fraction of center-of-mass energy trapped in the BH, y(z), by x m = M^ lin /[sy 2 (z)]. TeV BHs may carry 
electric or color charge and angular momentum. Immediately after formation, they are expected to decay through 
loss of excess multipole moments (balding phase), gravitational + Hawking radiation (37| (evaporation phase) and 
final n-body decay or remnant production (Planck phase) . SM particles are emitted on the brane and can be detected 
[HI]. Since the balding phase is poorly understood, simulations neglect the energy loss in this phase. The description 
of the evaporation phase is also approximated; since emissivities of rotating BHs are not known for all fields, BH 
generators use greybody factors for non-rotating BHs [H, |4(j. In CATFISH, the total decay multiplicity is [l9[ 

(n + i)s E*c^r Ki 



where Ci are the degrees of freedom of species i, T-p i and are the relative emissivities of Ref. [40| . S is the initial 
entropy of the BH, and V s and 1Z S are spin-dependent normalization factors. A more detailed discussion of the 
evaporation and Planck phases of TeV BHs can be found in Refs. [l3| . 

Dilcpton production in BH events differ greatly from dilepton production in the MSSM. Unlike SUSY, there is no 
single process of dilepton production; dileptons are either produced by the BH directly or by the decay of heavier 
particles such as the Zq boson, tt pairs of a combination of the two. Therefore, the BH dilepton invariant mass does 
not show a sharp cut-off at high energy. 



IV. EVENT SIMULATIONS 



SUSY simulations are carried out using a combination of ISAJET and PYTHIA, with the former generating the mass 
spectrum. BH simulations are carried out using the CATFISH Monte Carlo generator. The setup for each simulation 
is summarized below: 

• SUSY: 

— The MSSM mass spectrum is generated with ISAJET (ver. 7.75); 

— The mass spectra in SLHA format is fed into PYTHIA; 

— All SUSY processes except SM Higgs production are simulated; 



1 Recently, a new BH generator (BlackMax) appeared in the literature |31|. BlackMax is supposed to include BH rotational effects. 
However, the gravitational loss for rotating BHs is artificially set to zero. Since the energy loss due to gravitons is enhanced by rotation 
and extra-dimensional effects, and cannot be neglected, BlackMax results for rotating BHs should not be trusted. 
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— Unstable SM particles and sparticles are hadronized or decayed with PYTHIA. 
• BHs: 

— The cross section for a BH event is calculated in the center-of-mass frame; 

— The initial BH mass is sampled from the differential cross section; 

— The BH is decayed through Hawking mechanism and final n-body event (or remnant); 

— Unstable quanta are hadronized or decayed with PYTHIA. 

The benchmark model for SUSY is LHC point A. The parameters for the BH benchmark model are fundamental 
Planck scale M* = 1 TeV, minimum BH mass M m i n = 2 TeV, classical-to-quantum threshold Qmin — 1 TeV, six 
extra-dimensions (n = 6) and two-body final decay (n p — 2). Particles produced in the initial-radiation phase are 
removed by imposing Pt cuts of 5 GeV and 15 GeV for leptons and photons+hadrons, respectively (l9| . 



V. EVENT ANALYSIS 

In this section, we first use event shape variables to discriminate SUSY and BHs models. We then complement 
these results by looking at isolated dilepton events. The salient features of this analysis is that BH events tend to be 
more spherical than SUSY events due to the spherical nature of the Hawking radiation. This is specially evident for 
high-mass BHs. The formation of a stable BH remnant at the end of the evaporation phase also helps to discriminate 
MSSM and BH events because of the large amount of energy which is carried away by the remnant. Isolated dilepton 
events provide a further powerful means to distinguish the two models. This is due to the fact that leptons are rarely 
emitted by BHs (the hadron-to-lepton ratio is approximately 5:1) and are uncorrelated; they can be emitted at any 
angle w.r.t. beam axis, whereas SUSY dileptons typically originate from a single decay chain. 



A. ]f T and event shape variables 

Figure [3] shows visible energy, missing transverse momentum, and transverse momentum of leptons and 
hadrons+photons for 10,000 MSSM and BH benchmark events. Even in the absence of a BH remnant, the amount 
of visible energy and f T is comparable for the two scenarios. This is due to the presence of invisible channels in both 
models: the LSP for SUSY and neutrinos+gravitons for the BH. The flavor of the decay products is a better discrim- 
inator. MSSM interactions do not produce leptons with energy above the TeV since isolated leptons are produced 
by the decay of sparticles with typical energy of less than a few hundred GeV. On the contrary, quanta produced in 
the BH decay are characterized by an average energy E ~ M/N, where the multiplicity N is less than 10 for typical 
BHs at the LHC. Since Hawking evaporation does not distinguish leptons from hadrons, hard leptons with energy 
up to several TeV are likely to be produced during the BH decay. This suggests that isolated leptons may provide a 
powerful means to discriminate the two models. This is indeed the case, as we shall see in the next section. 

The four plots in the left panel of Fig. [4]show how variations in the BH Planck phase affect the observables of Fig. [3] 
The plots compare two- and four-body decays to the formation of a BH remnant (n p — 0). By the time the remnant 
has formed, the BH is expected to have shed electric and color charges. (See, however, Ref. 41] for an alternative 
scenario.) If this is the case, the BH remnant is undetectable and a source of missing transverse momentum in 
addition to neutrinos and gravitons which are emitted during the Hawking evaporation phase. This leads to a larger 
difference in f T between the MSSM and BH models. The visible transverse momentum in hadrons+photons is sensibly 
reduced in the presence of a BH remnant; the latter carries away energy which otherwise would have been emitted in 
visible channels (mostly hadrons) during the BH decay phase. It is interesting to note that the amount of transverse 
momentum in the leptonic channel is essentially unaffected by the presence of a BH remnant. This is due to the fact 
that leptons are rarer than hadrons in the BH decay phase; variations in the energy distribution of the leptonic channel 
are thus suppressed compared to the hadronic channel. Changes in the number of final Planckian hard quanta do not 
produce significant differences in the distributions; more quanta of lower energy behave statistically like less quanta 
with higher energy. Provided that the BH decays at the end of the Hawking phase, it is thus safe to set the number of 
Planckian quanta to n p — 2 or n p — 4, although BHs may decay in different numbers of particles on a event-to-event 
basis. Variations in the classical-to-quantum threshold Q m in are also not expected to cause significant differences in 
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FIG. 3: Comparison of 10,000 SUSY and BH benchmark events at the LHC. Visible energy and missing transverse momentum 
Ij> T (top panels) are comparable due to the presence of invisible channels in both models. Leptons with large transverse 
momentum provide instead an effective discriminator (bottom right panel). 



the energy /momentum distributions. A higher threshold increases the emission in the Planck phase while decreasing 
Hawking radiation. Since these phases differs only in relative greybody factors, the effect is too small to be detected. 

The four plots in the right panel of Fig. |4]show visible energy, ]f Tl and visible transverse momenta of leptons and 
hadrons+photons for different values of the fundamental Planck scale. Higher values of M* lead to more massive BHs, 
i.e. higher multiplicity and more energetic quanta. This causes a significant increase in missing and visible momenta. 
If the value of the fundamental Planck scale happens to be large, BHs are likely to be found and easily distinguished 
from SUSY through detection of highly-energetic isolated leptons and hadronic jets. Missing transverse momentum 
of several TeV would also be observed. 

Event shape variables such as sphericity and 2 nd Fox- Wolfram moment can be used to complement the above 
analysis. BH events arc more spherical because of the nature of Hawking radiation and the "jetty" nature of SUSY 
decays. Formation of a BH remnant and high values of the fundamental scale lead to significant higher sphericity 
than SUSY (top panels of Fig. [5]). The 2 nd Fox- Wolfram moment (bottom panels of Fig. [5]) is stable versus changes in 
the BH Planck phase and provides a good MSSM/BH discriminator. BH models with higher M* can be differentiated 
more easily from the MSSM. 

Similar conclusions can be reached by looking at jet masses and number of jets. The MSSM generates more and 
lighter jets than the BH model due to copious production of quarks (Fig. [5]). The difference is again specially significant 
for high values of M* and in the presence of BH remnants. Absence of sub-Q m i„ hard jets could provide strong evidence 
for BH remnant production. (See the suppression of heavy jets below the classical-to-quantum threshold Q m in = 2 
TeV in the top leftmost panel of Fig. [6]) 



B. Event analysis using high- Pt dileptons 



The use of isolated dileptons as SUSY signature has been extensively discussed in the literature [H, [42j ■ Although 
their production is not as high as colored particles, high-energy isolated leptons provide a cleaner environment by 
allowing the removal of the QCD background. Moreover, since most of BHs produced at the LHC are expected to be 
very light, multi-particle analysis may not provide the most effective discriminators (l9l l43j . The study of leptonic 
final states alleviates this problem. 
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FIG. 4: Distribution of visible energy, ]p T and transverse momenta of leptons and hadrons+photons. SUSY plots are shown 
as pink open squares. The four plots in the left panel show the effect of different decay modes in the Planck phase of ten- 
dimensional BHs: remnant formation (n p = 0, filled black triangles), two-body decay (n v = 2, filled red circles) and four-body 
decay (n p = 4, filled green squares). The fundamental Planck scale is M t = 1 TeV. The four plots in the right panel show the 
effect of varying the fundamental Planck scale: M t = f TeV (filled black triangles) , M* = 2 TeV (filled red circles) and Af* = 3 
TeV (filled green squares) . The ten-dimensional BHs decay in two hard quanta at the end of the evaporation phase. 



The dominant MSSM interaction for opposite-sign, same-flavor (OSSF) dileptons at LHC point A is [27| 
with a branching ratio of 27%. The maximum dilcpton invariant mass for this interaction is 



Mmax 




~\ 1/2 



100 GeV . 



The background for this process is due to SM decays of W, Z bosons and top quarks. This background can be 
removed by applying suitable cuts on transverse momentum and sphericity of the leptons [22| : 

• P T/ >15 GeV, 1*7x1 < 2.5; 

• Isolation cut, J^i Pti < 7 GeV in a cone of R = 0.2, 



where Pti is the transverse momentum of the leptons, R = \J Ar/ 2 + A0 2 



V 



■ ln[tan(0/2)] is the pseudorapidity, 



and 4> and are the azimuthal and polar angles of the lepton w.r.t. beam axis, respectively. 

Isolated leptons in BH decays come directly from the BH itself, from the decay of Zq bosons and top quarks, or 
from a combination of the two. Since the branching ratio of Zq into leptons is small, T(l + l^)/r tot ~ 0.034 [44], and 
the decay of top quarks into leptons is rare [45[, production of OSSF dileptons is less frequent in the BH model than 
in the MSSM. our analysis shows that an OSSF dilepton event occurs approximately every 100 BH and 20 SUSY 
events, with a ~ 1:5 ratio of BH-to-SUSY dilepton events at fixed luminosity. 

Figure [7] shows the dilepton invariant mass distribution for the MSSM (shaded pink plot) and the BH model with 
final two- and four-body decay (left and right panels, respectively). As expected, the SUSY distribution shows a sharp 
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FIG. 5: Sphericity (top panels) and 2 nd Fox- Wolfram moment (bottom panels) for 10,000 BH and MSSM events (pink open 
squares). The left panels show the effect of different Planckian decay modes: BH remnant (filled black triangles), two-body 
decay (filled red circles) and four-body decay (filled green squares). The fundamental scale is M* = 1 TeV and the number of 
extra dimensions is six. The right panels show the effect of different fundamental scales: Af* = 1 TeV (filled black triangles), 
2 TeV (filled red circles) and 3 TeV (filled green squares). The ten-dimensional BHs decay in two quanta at the end of the 
Hawking phase. 
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FIG. 6: Light and heavy jet masses (four leftmost panels) and number of jets (right panels) for 10,000 BH and MSSM events. 
Symbols are like in previous figures. 



edge at ~ 100 GeV [22[. The BH invariant mass distribution is characterized by a peak at ~ 90 GeV and a long tail 
up to energy of several TeV. The peak is due to dileptons events produced from the decay of Zq bosons, the dominant 
channel for OSSF dilepton production in BH. The tail is originated by uncorrelated lepton pairs emitted directly by 
the BH or in top quark decays. The leptons are hard and the reconstructed dilepton mass can have super- TeV values. 
The BH invariant mass for the two-body Planckian decay shows a second, smaller peak at ~ 1 TeV. This occurs 
because the BH at the end of the Hawking phase may decay in two OSSF leptons, leading to a reconstructed dilepton 
invariant mass equal to Q min . If the Planckian decay is a four-body process, the BH mass at the end of the Hawking 
evaporation is distributed among four quanta. This produces a lower, smoother reconstructed invariant mass. 
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FIG. 7: Invariant mass distribution (in GeV) for 1000 SUSY and BH OSSF dilepton events. The SUSY distribution (shaded 
pink histogram) shows the typical endpoint due to the presence of the LSP. The high- Pt tail of the BH distribution is originated 
by uncorrelated lepton pairs emitted during the Hawking evaporation phase. The final BH decay is in two-quanta (left panel) 
or four-quanta (right panel). 



The number of isolated, high- Pt leptons can also be used to complement the dilepton analysis (left panel of Fig.[S]). 
SUSY events are capable of producing up to five isolated leptons from the cascade decay of heavy sparticles. Events 
with X2X2 or X1X2 m ay produce four or three isolated leptons, respectively [46l |. On the contrary, events with three 
or more isolated leptons are very suppressed in BH decays at the LHC energy. Although multilcpton events are 
rare, there is very little background and they could be effectively used to distinguish the MSSM and the BH model. 
Other effective discriminators can be constructed by looking at dilepton events with same sign and/or opposite- flavor 
leptons. The "democratic" nature of the BH decay makes all dilepton events roughly equally probable, whereas the 
MSSM favors same-flavor dileptons. Presence of hard opposite-flavor leptons is a clear indication of BH decay (right 
panel of Fig. [8]). Our analysis shows that 73% of SUSY dilepton events are OSSF, compared to only 50% in the BH 
model. Conversely, opposite-flavor events are twice more frequent in the BH model (40%) compared to the MSSM 
(21%). 



VI. CONCLUSIONS 



We have discussed and compared the signatures of the MSSM and the TeV-BH model at the LHC. A thorough 
analysis of event-shape variables and dilepton events has shown that it is possible to distinguish the two models. 
BH events are characterized by higher sphericity than SUSY processes. If a BH remnant is formed at the end of 
the evaporation phase, missing Pt and heavy jet mass are effective signatures to discriminate BH formation from 
the MSSM. Although event-shape variables alone cannot unequivocally discriminate between SUSY and BHs, their 
knowledge may prove useful when combined with the analysis of the leptonic channel. Isolated dileptons could provide 
the "smoking gun" for detecting BHs at the LHC. The BH dilepton invariant mass shows a tail at high energy which 
is absent in the SM or MSSM. This analysis can be further strengthened by looking at the number and flavor of 
isolated leptons. 
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FIG. 8: Left Panel: Histogram of the number of events with high-Pr leptons for 10,000 MSSM (pink filled squares) and BH 
interactions (black open triangles) . The number of BH events with three isolated leptons is smaller than the number of SUSY 
events by a factor of ~ 20. The probability of producing BH events with four or more leptons is virtually zero. Right Panel: Pt 
scatter plot for ~ 1000 isolated opposite-flavor dilepton events for SUSY (pink filled squares) and BHs (black open triangles). 
BH leptons are harder than SUSY leptons and show a larger spread in Pt- 
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